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We propose a new method of the /3-7 constraint for quadrupole deformation in antisym- 
metrized molecular dynamics (AMD) to describe various cluster and shell-model structures 
in the ground and excited states of light nuclei. We apply this method to N = 6 isotones, 
(N ■ 10 Be, 12 C, 9 Li, and "B, and find various structures as functions of the deformation parame- 

S_i ' ters, /? and 7. In these nuclei, shell-model-like structures appear in the small ft region, while 

a : cluster structures develop well in the large j3 region where various geometric configurations 

■ of clusters are obtained depending on the 7 parameter. For 10 Be and 12 C, we superpose the 
basis AMD wave functions obtained by the /3-j constraint method to calculate energy spec- 

l/") ' tra, and prove the advantages of the present method of the two-dimensional p-j constraint 

in the framework of AMD. 

[ §1. Introduction 

In light nuclei, the cluster aspect is one of the essential features, as well as the 
shell-model aspect. Owing to the coexistence of these two natures, namely, cluster 
and shell-model features, various structures appear in stable and unstable nuclei. 
(T) , 12 C is one of the typical examples where the cluster and shell- model aspects 

coexist. The ground state of 12 C is known to have mainly a shell- model feature of the 
P3/2 subshell closed configuration, whereas the well-developed 3a-cluster structures 
fSJ ■ appear in excited states. In the theoretical works on the 3a-cluster structures J3~tSD 

various configurations of the 3a-cluster structures were suggested in the excited 
ON ' states above the 3a threshold energy, for example, the a condensation of weakly 

■ interacting three a clusters in the state and the equilateral-triangular structure 
I of three a clusters in the state. Moreover, a linear-chainlike (or an obtuse-angle- 
triangular) structure of three a clusters in the state was suggested. 

Cluster structures have also been found in light neutron-rich nuclei such as Be 
isotopes. In 10 Be, the low-lying states are understood in a molecular 2a + 2n pic- 
turep^'DS where two a cores are formed and two excess neutrons occupy molecular 
orbitals around the 2a. In terms of a simple shell model, 10 Be is an N = 6 nu- 
cleus, and therefore, the p 3 / 2 subshell closure effect is also important, as well as the 
2a + 2n cluster feature at least in the ground state. This means that the cluster- 
shell competition is essential in unstable nuclei as well as stable nuclei, as argued in 
Ref.ISD. 

For theoretical investigations of such nuclei, it is necessary to describe the coex- 
istence of shell and cluster features systematically. However, many theoretical frame- 
works still have deficiencies in describing both the shell-model and cluster structures. 
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In fact, in the case of C, shell models can be used to describe low-lying shell-model 
states but they usually fail to describe high-lying 3a-cluster states. On the other 
hand, conventional cluster models are suitable for studying the 3a-cluster states, but 
it is not easy to reproduce well the detailed properties of low-lying shell-model states 
because a cluster breaking is not incorporated in the cluster models. 

A method of antisymmetrized molecular dynamics (AMD^P^'Q® is one of the 
frameworks useful for overcoming this problem. It was applied to 12 C and suc- 
ceeded to describe the shell and cluster features due to the flexibility of its wave 
functions.®'® Moreover, in the study of fermionic molecular dynamics (FMD), in 
which model wave functions are similar to those of AMD, the coexistence of shell 
and cluster features in 12 C was described successfully.® 

The AMD method has also been applied to various stable and unstable nuclei, 
and it has been proved to be one of the powerful approaches of describing vari- 
ous structures such as cluster structures and shell-model structures P®~^ There 
are some versions of the AMD, for example, the variation after parity and total- 
angular- momentum projections (VAPs),® the variation with the constraint on the 
quadrupole deformation f3 (/5 constraint AMD)J^~^ or the constraint on the clus- 
ter distances (ci-constraint AMD)P§) In principle, a basis AMD wave function is 
given by a Slater determinant of Gaussian wave packets, and excited states are de- 
scribed by superposition of Slater determinants. In practical calculations of excited 
states of light nuclei, it is important to prepare efficiently various cluster configura- 
tions including 2-body and 3-body clusterings as basis wave functions in the AMD 
framework. Moreover, in the study of unstable nuclei, further flexible model wave 
functions such as 2-body or 3-body cluster structures with surrounding valence nu- 
cleons will be required to describe possible exotic cluster structures in excited states. 

To study a variety of cluster structures and the coexistence of cluster and shell 
features in light unstable nuclei, we propose an extended method of constraint AMD 
to describe various cluster and shell structures. That is the two-dimensional con- 
straint with respect to the quadrupole deformation parameters, j3 and 7, which is 
expected to be efficient for preparing basis wave functions with various cluster config- 
urations. We call this method /J-7 constraint AMD. We expect shell-model structures 
to appear in the small (3 region, whereas developed 2-body or 3-body cluster struc- 
tures can be obtained for large /3. In the large /3 region, various configurations of 
cluster structures may appear depending on (3 and 7. 

The /3-7 constraint AMD may also be useful in the study of triaxial deformations. 
On the other hand, in the Hartree-Fock-Bogolyubov (HFB) calculations, the j3- 
7 constraint was adopted, for example, in Ref. [2Tj) . and the superposition of /3-7 
constraint wave functions has been performed recently by Bender and HeenenJ^S 
It was found that the triaxiality is important to reproduce the experimental data 
of 24 Mg in the HFB calculations. However, works on triaxial calculations with the 
superposition are limited, and it is still a challenging problem. Moreover, such mean- 
field approaches are not necessarily suitable for describing cluster structures. For the 
study of cluster features, it is important to apply the constraint to a framework 
that can describe cluster structures. 

In this paper, we applied the /3-7 constraint AMD to N = 6 isotones, 10 Be, 12 C, 
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9 Li, and n B to check the applicability of this method. We analyze the results and 
confirm that various structures appear as functions of the deformation parameters, 
j3 and 7, in the present framework. In particular, we focus on the coexistence of 
shell and cluster features. For 10 Be and 12 C, we also calculate the energy spectra of 
excited states by the superposition of the obtained basis wave functions and compare 
the results with the experimental data. We show that the /3-7 constraint AMD is 
useful for reproducing the energy spectra. A role of the 7 degree of freedom is also 
discussed. 

The content of this paper is as follows. In £j21 we explain the framework of the 
/3-7 constraint AMD. The calculated results are shown in Sj3j In SjH we discuss the 
effect of the triaxial deformation parameter 7. Finally, in SjSJ a summary and an 
outlook are given. 



§2. Framework of (3-"f constraint AMD 



We adopt a method of AMD with constraint. The frameworks of AMD and 
constraint AMD are described in detail, for example, in Refs. [T5|)- H7p . In this paper, 
we propose a two-dimensional constraint with respect to quadrupole deformation 
parameters. 

2.1. Wave function of AMD 

In the method of AMD, a basis wave function of an ^4-nucleon system \<P) is 
described by a Slater determinant of single-particle wave functions \<fi) as 



1 



det{|v?i),--- ,\va)} 



(2-1) 
>, spin part 



The i-th single-particle wave function \(fii) consists of the spatial part 
\Xi), and isospin part \n) as 

\<Pi) = |^>|Xi>|r»). (2-2) 
is given by a Gaussian wave packet whose center is located at 



The spatial part 
Zi/y/v as 



<r| 



, 3 

7T 



exp 



(2-3) 



where v is the width parameter and is taken to be a common value for all the single- 
particle Gaussian wave functions in the present work. The spin orientation is given 
by the parameter £j, while the isospin part \t{) is fixed to be up (proton) or down 
(neutron), 



\Xi) = Citl t ) +&4J I >, 
\n) = \p) or \n). 



(2-4) 
(2-5) 



In a basis wave function |#), {X} = {Z,£} = {Z±, £1, Z2, £2, • • ■ ,Za,£a} are com- 
plex variational parameters and they are determined by the energy optimization 
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using the frictional cooling method.®'^ As the variational wave function, we 
employ the parity-projected wave function 

|ff*)=P±|#) = i±£|*>. (2-6) 

Here, P is the parity transformation operator. We perform the variation for the 
parity-projected energy (<^ ± |i/|<P =l= ) / (<^ ± |<? ± ), where H is the Hamiltonian. After the 
variation, we project the obtained wave function onto the total-angular-momentum 
eigenstate. It means that the parity projection is performed before the variation, 
and the total-angular-momentum projection is carried after the variation. 

2.2. /3-7 constraint 

To describe various cluster and shell-model structures that may appear in the 
ground and excited states of light nuclei, we constrain the quadrupole deformation 
parameters, j3 and 7, and perform the energy variation with the constraints on the 
/3-7 plane. 

The deformation parameters, /3 and 7, are defined as 

, C0S75 v p^)-g>-^> , ( , 7) 
R 2 = f «* 2 ) + (y 2 ) + (z 2 )) ■ (2-9) 

Here, (O) represents the expectation value of the operator O for an intrinsic wave 
function \<P). x, y, and z are the inertia principal axes that are chosen as (y 2 ) < 
(x 2 ) < (z 2 ) and (xy) = (yz) = (zx) = 0. To satisfy the latter condition, we also 
impose the constraints (xy)/R 2 = (yz)/R 2 = (zx)/R 2 = 0. To obtain the energy 
minimum state under the constraint condition, we add the constraint potential V^ ons t 
to the total energy of the system in the energy variation. The constraint potential 
Konst is given as 



V'const =Vl [03 cos 7 - fa cos 7 ) 2 + (/3 sin 7 - /3 sin 70) 



( (xy)\ 2 , / {yz)\ 2 . ( 
V R 2 ) \R 2 ) \R 2 



(2-10) 



Here, r/i and 772 take sufficiently large values. After the variation with the constraint, 
we obtain the optimized wave functions ^^(/^o, 70)) for each set of parameters, 
(/3,7) = (A),7o). 

In the calculations of energy levels, we superpose the total-angular-momentum 
projected wave functions P'l 1K |^ ± (/3, 7)). Thus, the final wave function for the 
state is given by a linear combination of the basis wave functions as 

K i 
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The coefficients f n {f3i,~fi,K) are determined using the Hill- Wheeler equation 

5 ({^ItflO - Eni^l^)) = 0. (2-12) 

This means the superposition of multiconfigurations described by parity and total- 
angular-momentum projected AMD wave functions. In the limit of sufficient basis 
wave functions on the /3-j plane, it corresponds to the generator coordinate method 
(GCM) with the two-dimensional generator coordinates of the quadrupole deforma- 
tion parameters, (5 and 7. 

2.3. Hamiltonian and parameters 

The Hamiltonian H consists of the kinetic term and effective two-body interac- 
tions as 



Coulomb 
ij j 



(2-13) 



Kj 



i<j 



i<j 



where V£F DkTai , V^ pm " or :t and v^? oulomb are the central force, spin-orbit force, and 
Coulomb force. As the central force, we use the Volkov No. 2 interaction,^ 



V; 



central 



E 

k=i 



v k exp 



Lll 



(W + BP a - HP T - MP a P T 



(2-14) 



where v\ = -60.65 MeV, v 2 = 61.14 MeV, a\ = 1.80 fm, and a 2 = 1.01 fm. For the 
spin-orbit part, we used the spin-orbit term of the G3RS interaction,^ which is a 
two-range Gaussian with a projection operator P( 3 0) onto the triplet odd state, 



yspm 
ij 



orbit 



E 

k=l 



exp 



P( 3 0) 



1 + P a 1 + P T 



P( 3 0)L • S, 



(2-15) 
(2-16) 



where b± = 0.600 fm and 62 = 0.477 fm. 

We take interaction parameters, the Majorana exchange parameter M = 0.6 
(W = 0.4), the Bartlett exchange parameter B = 0.125, and the Heisenberg exchange 
parameter H = 0.125 in the central force, and u\ = —1600 MeV and u 2 = 1600 
MeV in the spin-orbit force. All these parameters are the same as those adopted in 
Refs. I12p and I25|) . except for a small modification in the strength of the spin-orbit 
force. They are adjusted to reproduce the a + a phase shift (M, W = 1 — M, a\, a 2 ), 
binding energy of the deuteron {B = H), and a + n phase shift (u\ = —u 2 , b±, b 2 ). 
We slightly modify only the strengths of the spin-orbit force from —u\ = u 2 = 2000 
MeV adopted in Ref.QSJ) to -u\ = u 2 = 1600 MeV to fit the 0^ energy of 12 C in the 
present calculations. By using these parameters, the deuteron binding energy in an 
exact calculation is 2.24 MeV, which is in agreement with the experimental value, 
2.22 MeV. When the deuteron wave function is approximated by one AMD wave 
function, the calculated binding energy is 0.67 MeV, which slightly underestimates 
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the experimental value, 2.2 MeV. The Coulomb force Vh is approximated by 
a sum of seven Gaussians. 

For the width parameter of single-particle Gaussian wave packets in Eq. (|2-3|) . 
we used the value v = 0.235 fm~ 2 , which is determined from a variational calculation 
for the ground state of 9 Be in Ref. [25]) . 



§3. Results 



We applied the (3-j constraint AMD to the N = 6 isotones, 10 Be, 12 C, 9 Li, 
and n B. In this section, we analyze the results of the /3-7 constraint and show 
that various structures including well-developed cluster ones appear as functions of 
deformation parameters, /3 and 7, in the present framework. In particular, we focus 
on the coexistence of shell and cluster features. For 10 Be and 12 C, we also show the 
energy spectra of the excited states calculated by the superposition of the basis wave 
functions and compare them with the experimental data. 

3.1. Energy surfaces 

First, we performed variational calculations with the (3-j constraint at a total of 
196 mesh points of the triangle lattice on the (3-j plane. We obtained energy surfaces 
as functions of (3 and 7. The calculated energy surfaces on the f3--y plane before and 
after the total-angular-momentum projection for 10 Be, 12 C, 9 Li, and n B are shown 
in Figs. [H El [31 and|H respectively. We also show the energy curves along the axial 
symmetric line, 7 = 0° and 7 = 60°, in Fig. [5j 

In Fig. [1] for 10 Be, the left panel shows the energy of the positive-parity states 
before the total-angular-momentum projection, and the right panel shows the re- 
sults for + states calculated using the total- angular-momentum projection after 
the variation. We call the former the positive-parity energy surface and the latter 
the + energy surface. In the positive-parity energy surface, the minimum and a 
local minimum exist in the prolate region along the 7 = 0° line. The minimum 
point is at (/3 cos 7, /3 sin 7) = (0.35,0.0) in the normal deformation region, while 
the local minimum exists around {(3 cos 7, (3 sin 7) = (0.9,0.0) in the large deforma- 
tion region. After the total-angular-momentum projection, the shape of the energy 
surface changes because the energy gain due to the projection tends to be large in 
the large (3 region and triaxial region compared with that in the spherical or axial 
symmetric regions. As a result, the minimum point of the + energy surface shifts 
to {(3 cos 7, (3 sin 7) = (0.55,0.09) in the finite 7 region. The local minimum point 
of the + energy surface also has a finite 7 value as (f3 cos 7, f3 sin 7) = (1.03,0.04). 
Around the local minimum, the energy surface is rather flat. As we show later, the 
largely deformed excited band appears from this region in the GCM energy spectra. 

The calculated energy surfaces of 12 C are shown in Fig. [21 The positive-parity 
energy surface and the + energy surface are displayed in the left and right pan- 
els, respectively. The minimum point of the positive-parity energy surface is at 
{(3 cos 7, (3 sin 7) = (0.13,0.22), which indicates an oblate deformation. In contrast 
to 10 Be, which energetically favors the prolate deformation for any /3, in the case 
of 12 C, the prolate deformation around ((3 cos 7, f3 sin 7) = (0.8,0.0) is unfavored. 
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Fig. 1. Energy surfaces of Be on the j3-*f plane. The left shows the energy for the positive-parity 
states before the total-angular-momentum projection and the right shows that for the + states 
after the total-angular-momentum projection. 
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Fig. 2. Energy surfaces of 12 C on the /3-j plane. The left shows the energy for the positive-parity 
states and the right shows that for the + states after the total-angular-momentum projection. 



There is a peak with 7 MeV height at this point, and a valley goes from the 
oblate energy minimum to the prolate region keeping away from this peak. Af- 
ter the total-angular-momentum projection, the minimum point of the + energy 
surface becomes (/3 cos 7, /3 sin 7) = (0.35,0.17). This indicates that the deforma- 
tion of the energy minimum state changes from the oblate shape to the triaxial 
shape after the total-angular-momentum projection. There is also the flat region 
around (/3 cos 7, /3 sin 7) = (1.0,0.1) in the + energy surface. In the region around 
(/3 cos 7, [3 sin 7) = (0.6,0.0), the + energy surface has a large discontinuity because 
the intrinsic states in this region contain very small + components. 

The energy surfaces of 9 Li are shown in Fig. [3l The left and right panels show the 
negative-parity energy surface before the total-angular-momentum projection and 
the 3/2~ energy surface after the total-angular-momentum projection, respectively. 
In the negative-parity energy surface, it is found that the energy minimum state has 
a small deformation. A local minimum exists around (/3 cos 7, (3 sin 7) = (1.0,0.0) in 
the large prolate region. After the total-angular-momentum projection, the minimum 
point of the 3/2~ energy surface becomes (/3 cos 7, j3 sin 7) = (0.38,0.13). In the 
region near (/3 cos 7, /3 sin 7) = (1.03,0.04), the energy surface is rather flat, and 
there exists a local minimum. It is interesting that the behavior of the 3/2~ energy 
surface of 9 Li is qualitatively similar to that of the + energy surface of 10 Be. 

The results of n B are shown in Fig. 21 The negative-parity energy surface in the 
left panel shows the energy minimum point at (/3 cos 7, /3 sin 7) = (0.13,0.13). After 
the total-angular-momentum projection, the minimum point of the 3/2~ energy 
surface becomes (f3 cos 7, /3 sin 7) = (0.33,0.13) as seen in the right panel of Fig. 2J 
In the large prolate region, a valley is found around (/3 cos 7, /3 sin 7) = (0.9,0.1) in 
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Fig. 3. Energy surfaces of Li on the /3-7 plane. The left shows the energy for the negative-parity 
states and the right shows that for the 3/2~ states after the total-angular-momentum projection. 
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Fig. 4. Energy surfaces of n B on the /3~7 plane. The left shows the energy for the negative-parity 
states and the right shows that for the 3/2~ states after the total-angular-momentum projection. 
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Fig. 5. Positive-parity energy curves on the symmetric axes, 7 = 0° and 60° axes. The energy of 
the prolate states on the 7 = 0° line is displayed in positive f3 and that of the oblate states on 
the 7 = 60° line is displayed in negative j3. 



a similar way to 12 C. 

3.2. Structures on the fi-j plane 

In this section, we discuss the intrinsic structures obtained using the con- 
straint AMD, while paying attention to the cluster aspect. 

We analyze the spatial configurations of the Gaussian centers {Zi, Z2, • • • , Za} 
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and proton and neutron density distributions of each intrinsic wave functions 7)). 
We also investigate the difference p n — p p between the neutron density p n and proton 
density p p for neutron-rich nuclei to observe valence neutron behaviors. Here, the 
density p is that integrated along the y-axis as 



The density distributions for 10 Be are illustrated in Fig. El Figure Eta) shows 
the density distributions for \<P((3,~/)) at ((3 cos 7, (3 sin 7) = (0.35,0.00), which is the 
energy minimum point of the positive-parity energy surface. In this state, two a 
clusters are formed as seen in the dumbbell shape of the proton density. On the 
other hand, excess neutrons are found to occupy the axial symmetric orbitals and 
regarded to be in the -^^-like molecular orbitals. In fact, the expectation value 
of the squared neutron spin (S 2 ) is 1.07 and is close to the value (S 2 ) = 1 for 
the (vr 3 /2) 2 configuration. Figure 0(b) shows the density distribution of the energy 
minimum state with {f3 cos 7, (3 sin 7) = (0.55,0.09) in the + energy surface. It is 
found that two a clusters develop further. In this state, the density distribution 
of excess neutrons is not axial symmetric but is extended to the x-axis direction 
resulting in the axial asymmetric. This means that excess neutrons occupy the 
p x -like orbitals to have the finite 7 value instead of the axial symmetric molecular 
orbitals. Furthermore, (5^) becomes small to be 0.47. This indicates that excess 
neutrons have the larger S n = component than in the Fig. Eta) state and is 
consistent with the dominant (p x ) 2 configuration. This indicates that, after the total- 
angular-momentum projection, the energy minimum state has the more developed 
2a cluster structure and more two-neutron correlations of a spin-zero pair than 
those before the projection. Figures Etc), (d), and (e) are the density distributions 
for typical structures with oblate, prolate, and triaxial deformations in the large 
(3 region. Interestingly, two a clusters develop well in all these states with the 
large deformations, but differences are found in the distributions of excess neutrons. 
In Fig. Etc) for the oblate state, excess neutrons distribute far from the 2a core 
indicating a dineutron cluster. In Fig. EKd) for the almost prolate state, excess 
neutrons distribute around one of the two a cluster causing the 6 He correlation. In 
terms of the molecular orbital model, this state is also regarded as the molecular 
structure with the excess neutrons in the cr^-like orbitals as discussed in Ref . |26|) . 
Figure Ete) for the triaxial state shows two a clusters and a dineutron cluster with 
a configuration other than that in the oblate state. When we pay attention to the 
neutron densities of these three states, it is found that the neutron structure changes 
from the triangular configuration in the oblate region (c) to the narrow chainlike 
structure in the prolate region (d) as 7 becomes small. Therefore, the structures of 
10 Be as functions of deformation parameters, (3 and 7, are understood as follows. 
As (3 increases on the 7 = 0° axis, a cluster structure with two a clusters and 
surrounding two neutrons develops, and the excess neutron orbitals change from 
the 7r 3 /2-like orbitals to the cii/2-hke orbitals. On the other hand, as 7 becomes 
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(d) (1.03, 0.04) 
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Fig. 6. Density distributions of 10 Be. The proton density p p , neutron density p n , difference be- 
tween the neutron and proton densities p„ — p v integrated along the y-axis are illustrated 
in the left three panels. The density distributions of the intrinsic wave functions at (a) 
(P cos 7, P sin 7) = (0.35,0.00), (b) (/3 cos 7, /3 sin 7) = (0.55,0.09), (c) (/3 cos 7, p sin 7) = 
(0.28,0.39), (d) (P cos 7, /J sin 7) = (1.03,0.04), and (e) (/? cos 7, sin 7) = (0.70,0.26) on the 
/3-7 plane are shown. For the states (a) and (b), the difference between neutron and proton 
densities p' n — p' p , which is integrated along the z-axis but not the y-axis, is illustrated in the 
right panels. The size of the box is 10 x 10 fm 2 . 



finite from 0° to 60°, two neutrons have more the S n = correlation. When the 
deformation is sufficiently large, a dineutron cluster develops and drops out of the 
2a core. 

Next, we show the density distributions of 12 C in Fig. [7J In all the obtained 
wave functions |^(/3,7)), the proton and neutron densities are almost the same 
with each other. Figure [7Ja) is the density distribution for the minimum point 
(/3 cos 7, (3 sin 7) = (0.13,0.22) in the positive-parity energy surface. In this wave 
function, the centers of the single-particle Gaussian wave packets gather around 
the origin and there is no spatially developed cluster structure. The expectation 
value of squared intrinsic spin of protons and that of neutrons are 0.45, which 
indicates a large component of non-3a configurations. In other words, this state 
is considered to be the shell-model- like structure with the dominant P3/2 subshell 
closed component. Figure E|b) is the density distribution for the minimum point 
(/3 cos 7, j3 sin 7) = (0.35,0.17) in the + energy surface. In this state, an a cluster 
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Fig. 7. Density distributions of C. The proton density p v and neutron density p n are illustrated in 
the left and right columns, respectively. The density distributions of the intrinsic wave functions 
at (a) (j3 cos 7, /3 sin 7) = (0.13,0.22), (b) (fi cos 7, /3 sin 7) = (0.35,0.17), (c) (/3 cos 7, /3 sin 7) = 
(0.33,0.48), (d) 09 cos 7, /3 sin 7) = (1.13,0.04), and (e) (/3 cos 7, /3 sin 7) = (0.65,0.26) on the 
/3-7 plane are shown. The size of the box is 10 x 10 fm 2 . 



begins to develop compared with the state (a) for the energy minimum before the 
total-angular-momentum projection. However, single-particle Gaussian wave packets 
still gather around the origin, and this state is regarded as the intermediate between 
the shell-model-like structure and the cluster structure. In the large deformation 
region, three a clusters develop well in 12 C. As is expected, various configurations 
of 3 a clusters appear, depending on the deformation parameters, f3 and 7. Fig- 
ures [3(c), (d), and (e) are typical density distributions for the oblate, prolate, and 
triaxial deformed states, respectively. It is found that the equilateral-triangular con- 
figuration, linear-chainlike structure, and obtuse-angle-triangular configuration arise 
in the oblate state (c), prolate state (d), and triaxial state (e), respectively. Thus, on 
the /3-7 plane for 12 C system, we obtained various structures such as the shell-model 
structure in the small (3 region and the cluster structure in the large (3 region. In 
particular, three a clusters develop as the deformation parameter (3 becomes large. 
Various spatial configurations of three a clusters such as the linear chain and the 
equilateral-triangular structures are obtained as a function of the triaxiality 7. 

Figure [5] shows the density distributions of 9 Li. Figures [5]^ a) and (b) correspond 
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Fig. 8. Density distributions of 9 Li. The proton density p p , neutron density p n , and differ- 
ence between the neutron and proton densities p n — p P are illustrated in the left, middle, 
and right columns, respectively. The density distributions of the intrinsic wave functions at 
(a) (/?cos7,/3sin 7 ) = (0.15,0.00), (b) (/3 cos 7, /3 sin 7) = (0.38,0.13), (c) (/3 cos 7, P sin 7) = 
(0.28,0.39), (d) (J3 cos 7, P sin 7) = (1.03,0.04), and (e) (P cos 7, P sin 7) = (0.70,0.26) on the 
/3-7 plane are shown. The size of the box is 10 x 10 fm 2 . 



to the energy minimum states before and after the total-angular-momentum pro- 
jection, respectively. Figures [H^c) , (d), and (e) are density distributions of typical 
structures in the oblate, prolate, and triaxial states with large deformations, respec- 
tively. The cluster features in 9 Li are a good analogy to those in 10 Be, and therefore, 
they can be understood in a similar way as the above discussion for 10 Be by replacing 
one a cluster in the 10 Be system to a triton cluster in 9 Li. Namely, as the defor- 
mation parameter (3 increases on the 7 = 0° axis, the structure changes from the 
spherical shell-model state to the developed cluster structure with the a + 1 cluster 
and surrounding two excess neutrons in the axial symmetric molecular orbitals. As 
the triaxiality 7 increases from 0° to 60°, a pair of two neutrons drops out of the 
a + t cluster and describes the two-neutron correlation. 

The density distributions of n B are illustrated in Fig. The shell and cluster 
features of n B have a good correspondence to those of 12 C. In the small deformation 
region, the neutron density has an almost spherical shape and is regarded as the 
dominant p3/ 2 -shell closed configuration. As the deformation parameter (3 increases, 
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Fig. 9. Density distributions of 11 B. The proton density p p , neutron density p n , and differ- 
ence between the neutron and proton densities p n — p p are illustrated in the left, middle, 
and right columns, respectively. The density distributions of the intrinsic wave functions at 
(a) (P cos 7, P sin 7) = (0.13,0.13), (b) (P cos 7, P sin 7) = (0.33,0.13), (c) (P cos 7, P sin 7) = 
(0.33,0.48), (d) (P cos 7, P sin 7) = (1.13,0.04), and (e) (P cos 7, P sin 7) = (0.65,0.26) on the 
/3-7 plane are shown. The size of the box is 10 x 10 fm 2 . 



the structure changes from a shell- model state to a developed cluster structure via the 
shell-cluster competition region. The cluster features in the large deformation region 
are understood with the 2a + t clustering in 11 B instead of the 3a cluster structure 
in 12 C. Similarly to the case of 12 C, as a function of 7, we obtained various types 
of spatial configurations of three-center clusters such as the equilateral-triangular 
structure in the oblate deformation and the linear-chain structure in the prolate 
deformation. 

We summarize the shell and cluster features of these N = 6 systems obtained 
by the present method of the /3-7 constraint AMD on the f3-j plane. In the small (3 
region, /3 < 0.2 ~ 0.3, shell-model-like structures appear. In particular, the neutron 
structure has the dominant P3/2-shell closed configurations with the spherical shape. 
In this region, the centers of the single-particle Gaussian wave packets gather around 
the origin to keep the deformation small and there are no spatially developed cluster 
structures. In the large f3 region, cluster structures develop well. Depending on the 7 
parameter, various cluster configurations appear. That is, in the prolate region near 
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Fig. 10. Energy levels of the positive-parity states in Be. Four columns on the left are the 
experimental data, four columns in the middle are the two-dimensional /3-7 GCM results, and 
three columns on the right are the axial symmetric GCM results. In the two-dimensional /3-7 
GCM results, we classified the states into four groups. The first column is the ground band, the 
second column is the K n = 2 + side band of the ground band, and the third column is the 
band. Other states are plotted in the fourth column. The details are described in the text. 



the 7 = 0° line, two-body cluster structures or linear-chainlike structures develop 
well as P becomes large. In the oblate region near the 7 = 60° line, three-body cluster 
structures with equilateral-triangular or isosceles-triangular configurations appear. 
In the triaxial region, various triangle configurations of the three-body clusters are 
obtained. 

3.3. GCM results 

We superposed the wave functions obtained using the ^-7 constraint AMD with 
the GCM to calculate the energy spectra of 10 Be and 12 C. 

First, we describe the results for 10 Be. The calculated binding energy (B.E.) 
is 59.2 MeV, while the experimental one is 64.98 MeV. In the results, we obtained 
many excited states above the ground state. We show the calculated energy levels 
of 10 Be in Fig. \W\ as well as the experimental levelsP^~® 

The calculated results reproduce well the experimental energy levels. We clas- 
sified the calculated states in four groups by analyzing the components of the basis 
wave functions and intrinsic structures. 

The calculated 0^ , 2f, and 4f states constitute the ground band. They have the 
large component of the basis wave function at ((3 cos 7, /3 sin 7) = (0.55,0.09). For 
instance, the overlap of the + state with this wave function is 87%. As discussed 
before, in this dominant basis wave function, two a clusters develop well and excess 
neutrons occupy p x -like orbitals. The Of state also has a large overlap with the wave 
function at (/3 cos 7, j3 sin 7) = (0.35,0.00) with 70%. As already mentioned earlier, 
the wave function at (/3cos7, /3shi7) = (0.35,0.00) shows the structure with the 
2a core and excess neutrons in the ^^-like molecular orbitals. Therefore, roughly 
speaking, the ground bands can be interpreted as the molecular orbital structure 
with the 2a core and excess neutrons. However, it should be pointed out that 
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Table I. Electromagnetic transition strengths B(E2) and radii of 10 Be. (A) Calculated B(E2) 
values in the unit of e 2 fm 4 . The values in the GCM column are the two-dimensional 
GCM results and the values in the AS GCM column are the axial symmetric GCM results. (B) 
Root-mean-square radii for mass distributions of 10 Be calculated by the two-dimensional /J-7 
GCM. The unit is fm. 



(A) B(E2) (B) Root- mean-square radii 



Transitions 


/3- 7 GCM 


AS GCM 


Experiment 


States 


/3- 7 GCM 


Experiment 


2+^0+ 


9.4 


7.8 


10.24 ±0.97 


of 


2.39 


2.30 ±0.02 


2+^0+ 


1.2 


0.1 


0.64 ±0.23 


0+ 


2.98 




2+^0+ 


0.7 


0.1 




o 3 + 


2.96 




2+ -> 2+ 

z 2 ^ z l 


4.2 


0.8 











the dominant basis wave function is not the axial symmetric one but the triaxially 
deformed state with the excess neutrons in the p x -like orbitals. This means that the 
two-neutron correlation is contained significantly in the ground band. 

The 2\ and 3^~ states construct the K n = 2 + side band built on the ground 
band. These states appear because of the triaxiality of the ground band. Namely, 
the K 11 = 2 + band arises because the excess neutrons occupy the p x -like orbitals 
due to the two-neutron correlation. This is consistent with the theoretical work by 
Itagaki et al. where the triaxiality in 10 Be was discussedPD 

The 0^", 2^, and 4^ states belong to the excited K n = + band with a large 
prolate deformation. In fact, they have a large overlap with the wave function at 
(/3 cos 7, (3 sin 7) = (1.03,0.04), for example, as 64% overlap in the O^ state. As we 
discussed Fig. E^d), this dominant wave function has the extremely large prolate 
deformation due to the developed two a clusters and excess neutrons surrounding 
the 2a core. 

We also obtained other excited states below 14 MeV plotted in the rightmost 
column in Fig. [TOl The 0% state around 12 MeV is dominated by the basis wave 
function with developed two a and dineutron clusters. This might relate with the 
2a±dineutron condensate state suggested by Itagaki et 

The calculated E2 transition strengths are listed in Table [J (A) . The calculated 
strengths, B(E2, 2^ Of) = 9.4 e 2 fm 4 and B(E2, 2f -> Oj) = 1.2 e 2 fm 4 , agree well 
with the experimental data, B{E2, 2f -> Of) = 10.24 ± 0.97 e 2 fm 4 and B(E2, 2f -)• 
0+) = 0.64±0.23 e 2 fm 4 p3 respectively. The interband transition strengths between 
the ground band and its K w = 2 + side band are B{E2,2\ — > 0±) = 0.7 e 2 fm 4 and 
B(E2,22 — > 2^") = 4.2 e 2 fm 4 . The origin of these interband transition strengths 
is the recoil effect of the valence neutrons related to two a clusters as discussed in 
Ref.S]J. 

The calculated root-mean-square radii for mass distributions are listed in Ta- 
ble[I] (B). The calculated radius of the 0^ state, 2.39 fm, is close to the experimental 
one, 2.30 ± 0.02 fmP 1 

The present results of 10 Be are quantitatively similar to the earlier work studied 
with the molecular orbital modeP^ 1 and also to that with AMDp3 It indicates that 
the /?-7 constraint AMD is a useful method of preparing proper basis wave functions 
that are sufficient to describe the ground and excited states of 10 Be. 
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Fig. 11. Energy levels of the positive-parity states in C. Two columns on the left are the experi- 
mental data, three columns in the middle are the two-dimensional /3-y GCM results, and three 
columns on the right are the axial symmetric GCM results. In the two-dimensional GCM 
results, we classified the states in three groups. The first column is the ground band and the 
second column shows the spectra of the states, which have overlap with the basis AMD wave 
function in the broad 7 area of the large /3 region. The third column is the Ojj~ states. The 
details are described in the text. 



To check the advantages of the present method of the two-dimensional con- 
straint and to discuss the importance of the triaxiality, we also calculated the one- 
dimensional axial symmetric GCM in which the basis is limited to axial symmetric 
ones and compared the results with the two-dimensional GCM ones. This is a GCM 
calculation performed in a subset of two-dimensional /J-7 GCM. The number of su- 
perposed basis is only 40 on the 7 = 0° and 7 = 60° axes, whereas 196 basis are 
superposed in the two-dimensional f3-j GCM calculation. The calculated energy 
levels are shown in Fig. ITTJl The calculated B.E. is 58.6 MeV, which is just 0.6 MeV 
smaller than the /3-j GCM result. We found that axial symmetric basis wave func- 
tions are appropriate to reproduce the ground state. On the other hand, the axial 
symmetric GCM fails to reproduce some excited states. For example, the 2\ and 3^ 
states are 5.5 MeV higher than the /J-7 GCM result. Furthermore, the agreement of 
the calculated B(E2,2^ — > Of) with the experimental result is worse than the /3-7 
GCM result, as shown in Table HJ These failures are caused by the insufficiency of 
the GCM basis, which miss the triaxial ones. As a result of the lack of triaxial wave 
functions, interband transition strengths, B(E2^2^ — >• 0^) and B(E2^2^ — > 2^), 
become small. 

Next, we discuss the results of 12 C. The calculated B.E. is 92.5 MeV (The ex- 
perimental B.E. is 92.16 MeV). In Fig. [TT1 we give the calculated energy levels of 
12 C as well as the experimental data.® For the theoretical results, we classified 
the states in three groups by analyzing overlaps with the basis wave functions. In 
Fig. H21 the squared overlaps | (0 + |Pq( ) =0 |^ + (/3, 7))| 2 of the GCM wave functions for 
the + states with the basis AMD wave functions at each point (/3, 7) are shown. 

The 0^, and states constitute the ground band. They have a large overlap 
with the basis wave function at (/3 cos 7, (3 sin 7) = (0.35,0.17), for instance, as the 
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Fig. 12. Overlaps of the GCM wave functions for the + states of 12 C with each basis AMD wave 
functions obtained with the /3~7 constraint. The overlap with a certain (/?, 7) point is shown in 
the contour map. The solid lines show the j3 constraint paths, which indicate the optimum 7 of 
the energy minimum position for a given /? value. 



90% overlap in the Of state. As discussed earlier, this dominant basis wave function 
has no spatial development of three a clusters but the shell-model-like structure 
(Fig. [T](b)) . It is interesting that the dominant wave function shows the triaxiality 
in the intrinsic deformation. 

The calculated states displayed in the second column of the right side in Fig. [TT1 
have large overlaps with the basis AMD wave functions in the broad 7 area of 
the large /3 region. For example, as shown in Fig. [T2J the components of the 
state do not concentrate on some specific basis wave function but distribute widely 
into various basis wave functions. In this area of the large /3 region, the basis wave 
functions describe various configurations of the developed three a clusters depending 
on (3 and 7. It means that the Ojj" state is described using the linear combination 
of various 3a-cluster configurations. This is analogous to the earlier works with 
the 3a GCM,® FMDP and AMD calculations,® and is also consistent with the 
suggestion of the a condensate state.®' This indicates that the present method of 
/3-7 constraint is effective for describing the three-body cluster feature of the excited 
states described by a superposition of multiconfigurations. In particular, the degree 
of freedom of the triaxiality 7 is essential to describe this feature. A more detailed 
discussion will be given in the next section. 

The state has the largest overlap with the wave function at (/3 cos 7, ft sin 7) = 
(1.13,0.04) with 54% overlap. As shown in Fig. [T^d), this basis wave function has 
three a clusters with the linear-chain-like configuration. This result is consistent with 
the prediction of the 0^ state in the 3a model, FMD,® and AMD calculations P 

The calculated E2 transition strengths are listed in Table [TT] (A) . The calculated 
strengths, B(E2, 2f -)• 0^) = 6.0 e 2 fm 4 and B(E2, 2f -)• 0%) = 1.9 e 2 fm 4 , agree well 
with the experimental results, B(E2,2f 0^) = 7.59±0.42 e 2 fm 4 and B(E2,2f 
0^) = 2.6 ± 0.4 e 2 fm 4 fSJ respectively. 
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Table II. Electromagnetic transition strengths B(E2) and radii of 12 C. (A) Calculated B(E2) values 
in the unit of e 2 fm 4 . The values in the fi-q GCM column are the two-dimensional /3-~f GCM 
results and the values in the AS GCM column are the axial symmetric GCM results. (B) Root- 
mean-square radii for mass distributions of 12 C calculated by the two-dimensional /3-r/ GCM. 
The unit is fm. 
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The calculated root-mean-square radii for mass distributions are listed in Ta- 
ble [II] (B). The calculated value of the Of state, 2.31 fm, is close to the experimental 
value, 2.35 ± 0.02 fm or 2.31 ± 0.02 fmP^ 

The present calculations describe various structures in 12 C. The shell- model 
structure appears in the low-lying states, while 3a-cluster structures are found in 
the excited states. The shell and cluster coexistence obtained in the present results 
is consistent with the works with the FME0 and AMD calculations.® The success 
of the present calculations is due to the fact that the constraint AMD is effective 
for obtaining the basis wave functions for the ground and excited states of 12 C. 

Let us discuss the energies calculated by the axial symmetric GCM. The energy 
levels obtained by the axial symmetric GCM calculations are shown in Fig. [TTJ The 
calculated B.E. is 91.6 MeV, which is 0.9 MeV smaller than the f3-j GCM result. 
Since the difference from the /3-7 GCM result is small, the axial symmetric calculation 
is sufficient to describe the ground state, as in the case of 10 Be. On the other 
hand, the level structure of excited states changes considerably. For example, the 
0^ and 0% states are 2.6 and 3.4 MeV higher than the ,$-7 GCM result, respectively. 
Additionally, the B(E2, 2% — > 0^") value becomes small as shown in Table HT1 Clearly, 
the triaxial basis has an essential role in the description of excited states. 

§4. Discussion 

In calculations of deformed systems with the quadrupole deformation constraint, 
one often adopts only one-dimensional /? constraint instead of the two-dimensional 
/3-7 constraint. As mentioned earlier, various structures appear on the /3-7 plane 
and they are important basis wave functions, especially for the excited states. In 
this section, we discuss the advantages of the present method of the two-dimensional 
constraint and the importance of the triaxiality described by the 7 degree of freedom. 

The simplest assumption for deformed systems is the axial symmetry as carried 
out in deformed mean-field calculations. In such axial symmetric calculations, the 
deformation parameter f3 is the constraint parameter. This corresponds to the 7 = 0° 
and 7 = 60° lines on the f3-j plane. In Fig. [5j we have shown the energy curves on 
the 7 = 0°,60° lines for 9 Li, 10 Be, 11 B, and 12 C. The energy curves for the axial 
symmetric deformation show the double-well structure around f3 = 0. However, as 
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we mentioned earlier, it does not mean the coexistence of oblate and prolate shapes 
but it corresponds to one genuine minimum on the ^-7 plane for the ground state. 
For example, in the case of 10 Be, two local minima at (3 = —0.2 and j3 = 0.35 
are observed in the axial symmetric energy curve, as shown in Fig. [5j However, 
only the (3 = 0.35 minimum is the genuine minimum and the f3 = —0.2 minimum 
comes from the projection of the energy pocket of the genuine minimum at 7 = 0° 
onto the 7 = 60° line. In the case of 12 C, the local minima in the large prolate 
deformation around (3 = 1 is also artificial. This is a part of the valley in the region 
near (/3 cos 7, (3 sin 7) = (1.0,0.1). That is, a local minimum on the symmetric axis 
is sometimes an artifact owing to the axial symmetry restriction. Therefore, one 
should discuss carefully the energy curve of the axial symmetric calculations. In 
other words, the /3-7 constraint method is helpful for understanding the detailed 
behavior of the energy surface. 

We mention again the one-dimensional axial symmetric GCM. Needless to say, 
the basis wave functions with the triaxial deformations are missing in the axial sym- 
metric calculations, and therefore, some results for the excited states change consid- 
erably in comparison with the two-dimensional ^-7 GCM calculations. For example, 
the K w = 2 + side band of 10 Be cannot be described with the axial symmetric GCM 
calculations as discussed earlier. 

Another method of the one-dimensional constraint is the (3 constraint method 
without the axial symmetry assumption, which is usually applied to AMD calcula- 
tions. In this method, only the deformation parameter (3 is constrained, and the 
triaxiality 7 is a free parameter determined in the energy variation. Since 7 is au- 
tomatically optimized for each (3 constraint, it can take a nonzero value if a system 
favors a triaxial deformation. This method has been found to be useful for treating 
triaxial deformations of low-lying states fnjJHUJ >G23 Figure [TBI shows the paths of the 
j3 constraint method on the /3-7 plane. For reference, we also show the energy mini- 
mum point of the total-angular-momentum projected energy surface with the cross 
points in the figure. The energy curves along the paths obtained by the f3 constraint 
method are plotted as a function of (3 in Fig. [TH As seen in Fig. [TBI the (3 constraint 
path goes by the cross point for the energy minimum point of the total-angular- 
momentum projected energy surface but does not necessarily pass on the minimum. 
The difference in the minimum energy after the total-angular-momentum projection 
between the f3 constraint and /S-7 constraint is only 1 MeV in cases of 10 Be, 12 C, 
and 9 Li. Therefore, we can say that the f3 constraint method is reasonable, at least, 
in the description of the ground state properties. However, as we explain below, the 
basis wave functions obtained by the f3 constraint method are insufficient to describe 
some excited states of 10 Be and 12 C, and the two-dimensional ^-7 constraint is found 
to be essential. 

In the case of 10 Be, the prolate deformations are obtained mostly by the f3 
constraint, as seen in Fig. [TBI which shows that the (3 constraint path goes along the 
7 = 0° axis. In other words, the triaxially deformed structures in the finite 7 region 
cannot be obtained within the (3 constraint method. On the other hand, the 2^" and 
3f states in the K w = 2 + side band are obtained from the triaxially deformed basis 
wave functions. That is, the f3 constraint method is insufficient to reproduce the 
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Fig. 13. /3-7 values obtained by the {3 constraint AMD. The solid lines show the P constraint paths 
that indicate the optimum 7 ol the energy minimum position for a given j3 value. The energy 
surfaces before total-angular-momentum projection calculated by the fi-y constraint AMD are 
also shown with the contour map. Cross points stand for the energy minimum after the total- 
angular-momcntum projection obtained by the /3~7 constraint AMD. 
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Fig. 14. Energy curves calculated with /3 constraint AMD. 



K n = 2+ side band in 10 Be. 

In the case of 12 C, the two-dimensional constraint method is crucial for 
describing the 0^~ state as follows. We have shown the squared overlaps of the 
basis wave functions at each point (/3, 7) with the GCM wave functions for the + 
states in Fig. [T2J The Of and 0^~ states have large overlaps with the basis wave 
functions in the (/3 cos 7, j3 sin 7) ~ (0.35,0.17) and (/3 cos 7, j3 sin 7) = (1.13,0.04) 
regions, respectively. The (3 constraint line goes through these regions, and therefore, 
the 0^ and 0jj~ states may be reproduced with the basis wave functions on the /3 
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constraint line. However, the state has a large overlap with the basis wave 
functions in the broad 7 area of the large f3 region, which cannot be approached in 
the one-dimensional f3 constraint method. In this area, various configurations of the 
developed three a clusters appear depending on (3 and 7. As explained earlier, the 
linear combination of many 3a-cluster configurations is essential for describing the 
0^ state. These facts indicate that the /3 constraint method is insufficient and the 
/3-7 constraint method is necessary to describe the 0^" state of 12 C. 

§5. Summary and outlook 

To describe various cluster and shell structures in light nuclei, we proposed a 
new method of /3-7 constraint AMD. We applied this method to the N = 6 isotones, 
10 Be, 12 C, 9 Li, and n B. In the energy variation in the /3-7 constraint AMD, various 
structures appear depending on the two-dimensional constraint parameters, f3 and 
7. In the small {3 region, shell-model-like structures appear. In particular, the 
neutron structure has the dominant p 3 / 2 -shell closed configurations with the spherical 
shape. In the large f3 region, cluster structures develop well. As a function of the 
7 parameter, various cluster configurations appear. That is, in the prolate region 
along the 7 = 0° line, two-body cluster structures or linear-chainlike structures 
develop well as the deformation parameter f3 becomes large. In the oblate region 
along the 7 = 60° line, three-body cluster structures with equilateral-triangular or 
isosceles-triangular configurations appear. In the triaxial region, various triangle 
configurations of the three-body clusters are obtained. 

For 10 Be and 12 C, we superposed the basis AMD wave functions obtained by the 
/3-7 constraint method to calculate the energy spectra. We compared the present 
results with the experimental data. The results reproduce well the experimental 
spectra and are consistent with other theoretical studies. It was proved that the 
/3-7 constraint method is useful for describing the ground and excited states of these 
nuclei. We also mentioned the advantages of the two-dimensional /3-7 constraint 
method over the one-dimensional constraint method. It was confirmed that the 
degree of freedom of the triaxiality 7 in the two-dimensional /3~7 constraint has 
considerable importance to describe such excited states as the K w = 2 + side band 
of 10 Be and the 0^ state of 12 C. 

Moreover, in 9 Li and n B, cluster structures are expected to appear in the excited 
states from the similarity of the energy surface between 9 Li and 10 Be, and that 
between n B and 12 C. Further investigations with the GCM calculations using basis 
wave functions on the j3-j plane are required to discuss the energy spectra of 9 Li 
and U B. 

The present method proved to be very effective for describing various cluster and 
shell structures in the ground and excited states of light nuclei. We have applied the 
/3-7 constraint AMD to other nuclei such as 14 C, where shell- model- like features and 
cluster aspect are expected to coexist. We will report the results in a future paper. 
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